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FOREWORD 


NASA experience has indicated a need for uniform criteria for the design of space 
vehicles. Accordingly, criteria are being developed in the following areas of technology: 

Environment 

Structures 

Guidance and Control 
Chemical Propulsion. 

Individual components of this work will be issued as separate monographs as soon as 
they are completed. A list of all previously issued monographs in this series can be 
found on the last page of this document. 

These monographs are to be regarded as guides to design and not as NASA 
requirements, except as may be specified in formal project specifications. It is 
expected, however, that the criteria sections of these documents, revised as experience 
may indicate to be desirable, eventually will become uniform design requirements for 
NASA space vehicles. 

This monograph was prepared under the cognizance of the Langley Research Center. 
The Task Manager was A. L. Braslow. The authors were V. I. Weingarten and P. Seide 
of the University of Southern California and J. P. Peterson of NASA Langley Research 
Center. A number of other individuals assisted in developing the material and reviewing 
the drafts. In particular, the significant contributions made by E. H. Baker of North 
American Rockwell Corporation, C. D. Babcock, Jr., of California Institute of 
Technology, R. F. Crawford of Astro Research Corporation, J. B. Glassco of 
McDonnell Douglas Corporation, A. Kaplan of TRW Systems, M. H. Kural of Lockheed 
Missiles & Space Company, and J. Mayers of Stanford University, are hereby 
acknowledged. 

Comments concerning the technical content of these monographs will be welcomed by 
the National Aeronautics and Space Administration, Office of Advanced Research and 
Technology (Code RVA), Washington, D.C. 20546. 


August 1968 
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xy 


transverse shear-stiffness parameter for isotropic sandwich wall, 
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^XZ 1 

h — ^ (ti + t 2 ) 

bending stiffness per unit width of wall in x- and y-direction, 
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BUCKLING OF THIN-WALLED 
CIRCULAR CYLINDERS 


1. INTRODUCTION 


Structural components are said to be unstable under static loading when infinitesimal 
load increases or other small disturbances will induce the structure to change from one 
equilibrium configuration to another of a different character. For some structures and 
loadings the two configurations may differ slightly and large changes of shape develop 
gradually with an increase of load. In this case, the initial buckling load merely 
indicates when a change in the character of the deformations will occur. Generally, a 
more significant load is the ultimate load of the structure which may be reached either 
when the material fails plastically or when the structure collapses. In other instances 
the two equilibrium configurations differ greatly and the transition from one to the 
other is extremely rapid. The structure usually loses its capacity to sustain further load 
increases or it undergoes large deformations that render it unsafe for continued use. 

The primary design problem is the prevention of the buckling that leads to undesirable 
configurations — in particular, collapse. The magnitude of the critical static load ol a 
structure generally depends on its geometric proportions, the manner in which it is 
stiffened and in which it is supported, the bending and extensional stiffnesses of its 
various components, or other reasonably well-defined characteristics. For thin-walled 
shell structures, less certain characteristics, such as small deviations of the structure 
from its nominal unloaded shape, may also have quite important effects on the load at 
which buckling will occur. Other factors that affect buckling, such as cutouts, 
nonuniform stiffnesses, and variation of loading with time, are not considered in this 
monograph. 

This monograph indicates current practices for predicting buckling ol uniform 
stiffened and unstiffened circular cylindrical shells under various types of static 
loading, and suggests the procedures that yield estimates of static buckling loads 
considered to be conservative. The buckling of truncated conical shells and shells of 
double curvature will be treated in separate monographs. 


Estimation of design loads for buckling involves the use of the ultimate design factor. 
Considerations involved in selecting the numerical value of this factor will be presented 
in another monograph. 



2. STATE OF THE ART 


A relatively detailed understanding of the theory of buckling of circular cylindrical 
shells has been made possible by the use of electronic digital computers. Such 
understanding has been immeasurably aided by more rigorous formulations of the 
theory and by reliance on actual investigation, rather than on intuition. In particular, 
the large body of information available for flat plates is now known not to be strictly 
applicable to circular cylindrical shells, for which a separate body of information is 
being generated. 


Most of the available information on buckling of circular cylindrical shells is restricted 
to unstiffened shells of uniform thickness or to stiffened shells with uniform stiffness 
and properties, subjected to axisymmetric loading states which have certain simple 
longitudinal distributions, generally uniform. While problems involving surface loadings 
and thickness or stiffness variations that are axisymmetric but nonuniform 
longitudinally have been solved, detailed investigations of the effects of various 
parameters have not been made. Available information is quite inadequate for the 
analysis of loadings that are nonuniform circumferentially. Problems of this type can 
be treated by digital computer programs. A program for shells with uniform wall 
stiffnesses under axisymmetric loading is given in reference 1 and one for shells with 
axisymmetric geometric properties but asymmetric loadings in reference 2. 

Application of theory to the design of actual cylindrical shells is complicated by 
apparent discrepancies between theory and experiment. For shells in which 
longitudinal compression of the cylinder wall predominates, the discrepancies can be 
quite large. For shells in which shear or circumferential compression predominates, the 
discrepancies are generally less severe but still large enough to require experimental 
programs to establish design data. The causes of such discrepancies are generally 
understood. 

The primary source of error is the dependence of the buckling load of cylindrical shells 
on small deviations from the nominal circular cylindrical shape of the structure. 
Because the unloaded shape of test specimens usually has not been stringently 
controlled, most test results for nominally identical specimens have large scatter and 
fall below the theoretical values. 

Another source ot discrepancy is the dependence of buckling loads of cylindrical shells 
on edge values of longitudinal and circumferential displacements or forces. Because 
tangential edge conditions have not usually been precisely controlled in buckling tests, 
some of the scatter of test results can be attributed also to this source. Current 
methods of establishing design data tend to treat both initial imperfections and edge 
conditions as random effects. Results from all available tests are lumped without regard 
to specimen construction or methods of testing, and are analyzed to yield lower bound 
or statistical correction factors to be applied to simplified versions of the theoretical 
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results. This technique has proved satisfactory to date for design purposes. To put 
improved procedures to immediate use, however, the designer is advised to be alert to 
new developments in shell-stability analysis. The recommendations will be modified as 
more theoretical and test data become available. 


3. CRITERIA 
3.1 General 


Structural components consisting of thin, curved isotropic or composite sheet, with or 
without stiffening, shall be so designed that (1) buckling resulting in collapse of the 
structural component will not occur from the application of design loads, 
and (2) buckling deformation resulting from limit loads will not be so large as to 
impair the function of the structural components or nearby components, nor so large 
as to produce undesirable changes in loading. 


3.2 Guides for Compliance 

Design loads for buckling are considered to be any combination of ground or flight 
loads, including loads resulting from temperature changes, that cause compressive 
inplane stresses (multiplied by the ultimate design factor) and any load or load 
combination tending to alleviate buckling (not multiplied by the ultimate design 
factor). For example, external pressure loads or torsional loads should be increased by 
the design factor, but internal pressure loads should not. 

Suitably devised tests are required of representative structures under conditions 
simulating the design loads when minimum weight is a dominant factor or when 
cutouts, elastic end supports, or other special problems occur in the design. 


4. RECOMMENDED PRACTICES 
4.1 Scope 

Within the limitations imposed by the state of the art, acceptable procedures for the 
estimation of critical loads on circular cylindrical shells are described in this section. 
The important problems are indicated and the source of the procedures and their 
limitations are discussed. Where the recommended procedure is complex and is suitably 
defined in all its detail in a readily available reference, it is merely outlined. Where 
practicable, a summary of the procedure is given. 
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4.2 Isotropic Unstiffened Cylinders 


Unstiffened isotropic circular cylinders subjected to various conditions of loading are 
considered in this section. In the theoretical analysis of cylinders it is usually necessary 
to take account of deformations of the cylinder prior to buckling and of end 
conditions associated with displacements, both tangential and normal to the middle 
surface. The differences between rigorous solutions for various end support conditions, 
however, are obscured by the effects of initial imperfections. It is therefore customary 
to use simplified theoretical calculations modified to fit available test data. 

4.2.1 Axial Compression 

Buckling and collapse coincide for the isotropic circular cylinder in axial compression. 
From reference 3, an approximate buckling equation for supported cylinders under 
axial compression is 


where 


^x ^x 


7 T 2 D 


k x = m 2 (1 + )3 2 ) 2 + — 


2 7 2 
7 Z 


71 m 2 (1 + p ) 


(1) 


( 2 ) 


and the factor 7 2 has been inserted in the second term to correct the disparity between 
theory and experiment. Minimization of the equation with respect to m and j 3 results in 
the variation of k x with yZ shown in figure 1. For moderately long cylinders 
(yZ > , the buckling coefficient may be expressed approximately by the 

equation 


kx 



71 


(3) 


Equation (3) can be rewritten in a more familiar and useful form as the equation for 
axial stress: 


o 


X 


7^ _t_ 

Zhi - M 2 ) r 


= 0.67K-p(tor ju = 0.3) 


(4) 


4 



X 



2 468 2 468 2 468 2 468 2 468 

1 cf 1 1 10 10 2 10 3 10 4 


7Z 


Figure 1 

Buckling coefficients for simply supported isotropic circular 
cylinders subjected to axial compression 


The theoretical value of the factor y 2 in equation (2) is equal to unity in reference 3. A 
more rigorous analysis (ref. 4) indicates that lower theoretical values should be used for 
all types of support and that the shape of the curve is not correct in the low yZ range. 
On the basis of various experimental data (ref. 5), however, it is recommended that 
figure 1 be used with the factor y taken as 


7=1- 0.90 1(1- e'0) 


(5) 


where 


<t> = 77 /T for (7 < 1500) ( 6 ) 

16 v t t 

Equation (5) is shown graphically in figure 2 and provides a good lower bound for 
most test data. Figure 2 should be used with caution for cylinders with length-radius 
ratios greater than about 5 since the correlation has not been verified by experiment in 
this range. 
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Figure 2 

Correlation factors for isotropic circular cylinders subjected 
to axial compression 


When geometry and material properties are such that the computed buckling stresses 
are in the plastic range, the value E in equations (3) and (4) should be replaced by the 
value 77 E where 


V = 



tan 


(7) 


Equation (7) is an approximation of the plasticity factors given in references 6 and 7, 
and applies to moderately long cylinders. For extremely short cylinders (Z -► 0), the 
appropriate plasticity factor as given by reference 8 is: 


V = 


^tan 

E 


( 8 ) 


For cylinders with a length between those for which equations (7) and ( 8 ) apply, a 
plasticity factor is not available. Presumably, a linear interpolation with Z between the 
factors given by equations (7) and ( 8 ) would provide satisfactory results. Such a factor 
would be a function of cylinder geometry, as well as of the usual material stress-strain 
curve. 
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4.2.2 Bending 


Buckling and collapse coincide for isotropic, unpressurized circular cylinders in 
bending. The procedures given for isotropic cylinders in axial compression may be used 
also to obtain the critical maximum stress for isotropic cylinders in bending, except 
that a correlation factor based on bending tests should be used in place of the factor 
given by equation (5) for cylinders in axial compression. The correlation factor tor the 
cylinder in bending is taken on the basis of reference 9 as 

7=1- 0.731 (1 - e*) (9) 

where 



Equation (9) is presented graphically in figure 3. This equation should be used with 
caution for y > 1500 because experimental data are not available in this range. 
Although the theoretical critical stress for axial compression and that for bending are 
the same, the correlation factor for bending is greater than that for compression. This 



io 2 io 3 

r/t 


Figure 3 

Correlation factors for isotropic circular cylinders subjected 
to bending 
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is primarily because the buckling of a cylinder in compression can be triggered by any 
imperfection on the shell surface, whereas in bending, buckling is generally initiated in 
the region of the greatest compressive stress. 

4.2.3 External Pressure 

The term “lateral pressure” designates an external pressure which acts only on the 
curved walls of the cylinder and not on the ends. The load in the cylinder wall is given 
by 


Ny = Oyt = pr (11) 

The term “hydrostatic pressure” designates an external pressure which acts on both the 
curved walls and the ends of the cylinder. The cylinder wall loads are given by 

Ny = Oyt = pr (12a) 

N x = ° X t=y (12b) 

Except for unusually short cylinders, the critical pressures for the two types of loads 
are not significantly different. 

An approximate buckling equation for supported cylinders loaded by lateral pressure is 
given by reference 3 as 


where 



tt 2 D 

e 2 


k 


y 


prfi 2 

7T 2 D 


1 0 l2 



2~2 
7 Z 

(1 +/ 3 2 ) 2 


(13) 


(14) 


The buckling equation for cylinders loaded by hydrostatic pressure is obtained by 
replacing k y by k„ and the factor 0 2 before the bracketed expression in equation (14) 
by the factor (/3 + ). The factor y 2 has been inserted in the last term of 

equation (14) as a correction for the disparity between theory and experiment. 


The minimum values of k y for lateral pressure and k p for hydrostatic pressure are 
obtained by allowing the buckle aspect ratio /? to vary continuously and are shown in 
figure 4. The straight-line portion of the main curve (yZ > 100) is given by the 
equation (from ref. 3): 
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Buckling coefficient. 




Figure 4 

Buckling coefficients for simply supported isotropic circular 
cylinders subjected to external pressure 
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The family of curves at high values of 7 Z in figure 4, which is dependent upon the 
radius-thickness ratio of the cylinder, and which corresponds to buckling of the 
cylinder into an oval shape, is given by reference 10 as 



The theoretical value of 7 given in reference 3 is equal to unity. Other analyses (refs. 4 
and 11 ) indicate that restraint of the edge against longitudinal movement and/or 
rotation can increase the theoretical buckling pressure by as much as 50%. 

Experimental data for cylinders which buckle with more than two circumferential 
waves [ 7 Z < 11 . 8 (y) (1 — ju 2 )] show considerable scatter about the theoretical 
values given by equation (14) with 7 equal to unity. See reference 12. Some of this 
scatter is because the actual end restraint of test specimens is usually not considered in 
detail in the analysis of the test data. Scatter may also occur when isolated buckles 
appear in cylinders with large y or small y before a pressure is reached at which a 
buckle pattern appears around the entire shell circumference. The definition of 
buckling for these cases is a matter of individual judgment and may vary in different 
tests by different investigators. Because some tests fall as much as 25% below the 
theoretical results, a correlation factor of 

\fy = 0.75 (20) 

is recommended for use with equations (15), (16), and (17). 

For long cylinders that buckle into an oval shape, there is less discrepancy between 
theory and experiment (ref. 13), and a correlation factor of 

7 = 0.90 (21)_ 

is recommended for use with equations ( 1 8) and ( 1 9). 

For unusually short cylinders under lateral pressure ( 7 Z < 5), the plasticity factor for 
long, simply supported plates in axial compression may be used. It is obtained from 
reference 8 as 


Hsec /\ + l /l + 3 Etan\ 
F \2 2j4 4 E sec J 


( 22 ) 
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For 100 < 7 Z < 11.8 (1 - p 2 ), the approximate plasticity factor is obtained 

from equation ( 1 2 ) of reference 1 3 as 


V 



1 3 

— H 

4 4 



(23) 


and for yZ > 11.8 ({)" (1 -p 2 ), the approximate plasticity factor is obtained from 
equation (59) of reference 7 as 


V 


E 


sec 


E 




(24) 


No plasticity factor is available for the range 5 < 7 Z < 100; satisfactory results may 
however, be achieved by linear interpolation with the parameter Z between the values 
of 77 given by equations (22) and (23). 


Plasticity factors for the biaxial stress state of hydrostatic pressure are unavailable. For 
lack of better information, the plasticity factors given by equations (22) and (24) may 
be used. 


4.2.4 Torsion 


Buckling and collapse torques for unstiffened cylinders in torsion generally coincide. 
The theoretical buckling coefficient for cylinders in torsion can be obtained from 
figure 5, which is taken from reference 3. The straight-line portion of the curve is given 
by the equation 


"xy 


N X y£ 

tt 2 D 


= 0.85 ( 7 Z) 


(25) 


and applies for 50 < 7 Z < 78 (jY (1 -p 2 ). The correlation factor 7 has been 
inserted so that the theory may be adjusted to agree with experimental data. 
Equation (25) can be written 


r xy 


0.7477 4 H 


5 1 



(26) 
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> 

X 

-X 



o> 

c 



2 468 2 468 2 468 2 468 2 468 
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7Z 


Figure 5 

Buckling coefficients for simply supported isotropic circular 
cylinders subjected to torsion 


For yZ > 78 (-^j (1 — ju 2 ), the cylinder buckles with two circumferential waves. The 
buckling coefficient as given by equations (1 1) to (27) of reference 10 (p. 504) is: 


or 


^xy 


2/2 yZ 




yE 


'xy 


3 V7 (i — p 2 ) 4 


-(/ 


(27a) 


(27b) 


To approximate the lower limit of most data (ref. 3), the value 

2 

7 4 = 0.67 (28) 

is recommended for moderately long cylinders. 
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Plasticity may be taken into account by applying the plasticity factor from reference 
13 in the above equations: 


V = 


^sec 

E 


(29) 


The quantity E sec is obtained from a uniaxial stress-strain curve at a normal stress 
equal to twice the critical shear stress. Equation (29) applies to cylinders of all lengths. 


4.2.5 Combined Loads 


Load combinations usually encountered in practice are treated here. Generally, the 
recommended practice for handling combinations of two or more load conditions that 
may cause buckling is that the sum of the various critical load ratios is equal to unity. 


4. 2. 5.1 Circular Cylinders in Axial Compression and Bending 

The recommended interaction equation for combined compressive load and bending is 

R c + R b = 1 (30) 


The quantities R c and R^ are, respectively, the compressive and bending load or stress 
ratios. The denominators of the ratios are the allowable loads or stresses given by 
equations (4) to (6) for cylinders in axial compression and by equations (4) and (9) for 
cylinders in bending. 


4. 2. 5.2 Circular Cylinders in Axial Compression and External Pressure 

The recommended interaction equation for combined compressive load and 
hydrostatic pressure or lateral pressure is 

R c + R p = 1 (31) 

The quantities R c and Rp are, respectively, the compressive and hydrostatic- or 
lateral-pressure load or stress ratios. The denominators of the ratios are the allowable 
stresses given by equations (4) to (6) for cylinders in axial compression and by 
equations (16) and (20) for cylinders subjected to external pressure. 
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4.2. 5.3 Circular Cylinders in Axial Compression and Torsion 


For cylindrical shells under torsion and axial compression, the analytical interaction 
curve is a function of Z, and varies from a parabolic shape at low-Z values to a straight 
line at high-Z values. The scatter of experimental data suggests the use of a straight-line 
interaction formula. Therefore, the recommended interaction equation is 

R c + R t = 1 (32) 

The quantities R c and Rt are, respectively, the compressive and torsion load or stress 
ratios. The denominators of the ratios are the allowable stresses given by equations (4) 
to (6) for cylinders in axial compression and by equations (26) and (28) for cylinders 
in torsion. 


4.2. 5.4 Internally Pressurized Circular Cylinders in Axial Compression 


Buckling and collapse coincide for internally pressurized circular cylinders in axial 
compression as in the case of the unpressurized cylinder. Pressurization increases the 
buckling load in the following ways: 

1 . The total compressive load must be greater than the tensile pressurization 
load p7rr 2 before buckling can occur. 

2. The destabilizing effect of initial imperfections is reduced. The 
circumferential tensile stress induced by the pressurization inhibits the 
diamond buckling pattern, and, at sufficiently high pressurization, the 
cylinder buckles in the classical axisymmetric mode at approximately the 
classical buckling stress. 

Lower bound curves giving the increase in buckling load as a function of pressure, 
based on results for Mylar cylinders, are given in reference 14 for various 
radius-to-thickness ratios. Because these curves are unsubstantiated at present for other 
materials, the more conservative values given in reference 15 are recommended for 
design use. It is therefore recommended that the total load for buckling, unless 
substantiated by test, be obtained by the addition of the pressurization load p?rr 2 , the 
buckling load for the unpressurized cylinder [eqs. (4) and (5)] , and an increase in the 
buckling load caused by pressurization; that is 

Ppress = 2rrEt 2 /y 3(i 7 ~ + + Pm' 2 (33) 

where A7 is obtained from figure 6. 


For 11 - 0.3 
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Figure 6 

Increase in axial-compressive buckling-stress coefficient of 
cylinders due to internal pressure. 

4.2. 5. 5 Internally Pressurized Circular Cylinders in Bending 


For thin-walled cylinders subjected to bending and internal pressure, collapse loads are 
considerably higher than buckling loads (refs. 16 to 18), with the increase being 
substantially more than the tension stress induced by the pressurization. For example, 
for the true membrane cylinder (ref. 19, p. 229), the collapse load (M = p7rr 3 ) is twice 
the initial buckling load. The theoretical collapse load is, however, unattainable unless 
large undesirable deformations are present. It is therefore recommended that the 
collapse moment for pressurized cylinders be obtained by adding the moment-carrying 
capability of a pressurized membrane cylinder (taken for design purposes as 80% of the 
theoretical value), the collapse moment for the unpressurized cylinder [eqs. (4) and 
(9)], and an increase in the critical moment caused by pressurization. Then 
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press 


= TTlFt 2 
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(35) 


where A7 is obtained from figure 6. 

Forju = 0.3 

M press = (O.67 + A7) + 0.8 pwr 


(36) 


15 




42.5.6 Internally Pressurized Cylinders in Axial Compression and Bending 

For internally pressurized circular cylinders in combined axial compression and 
bending, equation (30) is recommended for use in combination with equations (33) 
or (34) and (35) or (36). 

4.3 Orthotropic Cylinders 

The term “orthotropic cylinders” covers a wide variety of cylinders. In its strictest 
sense, it denotes cylinders made of a single orthotropic material or of orthotropic 
layers. It also denotes types of stiffened cylinders for which the stiffener spacing is 
small enough for the cylinder to be approximated by a fictitious sheet whose 
orthotropic bending and extensional properties include those of the individual 
stiffening elements averaged out over representative widths or areas. Generally, the 
directions of the axes of orthotropy are taken to coincide with the longitudinal and 
circumferential directions of the cylinder. 

The behavior of the various types of orthotropic cylinders may be described by a single 
theory, the elements of which are equations of equilibrium for the buckled structure, 
relationships between force and moment resultants, and extensional and bending 
strains. For cylinders of a single orthotropic material, it is generally permissible to 
neglect coupling between force resultants and bending strains, and between moment 
resultants and extensional strains. The theory is then similar to that for isotropic 
cylinders. For stiffened cylinders or for cylinders having orthotropic layers, however, 
the neglect of the coupling terms can lead to serious errors. For example, the inclusion 
of coupling terms yields a significant difference in theoretical results for stiffened 
cylinder configurations having stiffeners on the inner surface or on the outer surface. 
The difference vanishes when coupling is omitted. 

Theoretical and experimental results for stiffened shells are generally in better 
agreement than those for unstiffened shells. The possibility of local buckling of the 
cylinder between stiffening elements should be checked. 


4.3.1 Axial Compression 


A buckling equation for stiffened orthotropic cylinders in compression (ref. 20) is 
given by: 
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Values of stiffnesses to be used for various types of construction are given in Section 
4.3.6. Prebuckling deformations are not taken into account in the derivation of the 
equation. The cylinder edges are assumed to be supported by rings that are rigid in 
their own plane but offer no resistance to rotation or bending out of their plane. The 
equation can be specialized for various types of cylinders which have been treated 
separately in the literature; for example, unstiffened or stiffened orthotropic cylinders 
with eccentricity effects neglected (refs. 21 and 22), and unstiffened or stiffened 
orthotropic cylinders with eccentricity effects taken into account (refs. 23 to 27). For 
ring-stiffened corrugated cylinders, a similar but not identical theory is given in 
references 28 and 29. For given cylinder and stiffener dimensions, the values of m and 
n to be used are those which minimize N x . 
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The unusually large number of parameters in equation (37) does not permit any 
definitive numerical results to be shown. For combinations of parameters 
representative of stiffened shells, calculations indicate that external stiffening, whether 
stringers or rings, or both, can be more effective than internal stiffening for axial 
compression. Generally, calculations neglecting stiffener eccentricity yield 
unconservative values of the buckling load of internally stiffened cylinders and 
conservative values of the buckling load of externally stiffened cylinders (ref. 23). An 
extensive investigation of the variation of the buckling load with various stiffener 
parameters is reported in reference 25. The limited experimental data (refs. 28 to 35) 
for stiffened shells are in reasonably good agreement with the theoretical results for the 
range of parameters investigated. 

On the basis of available data, it is recommended that the buckling loads of cylinders 
with closely spaced, moderately large stiffeners calculated from equation (37) be 
multiplied by a factor of 0.75. Correlation coefficients covering the transition from 
unstiffened cylinders to stiffened cylinders with closely spaced stiffeners have not been 
fully investigated. While theory and experiment (ref. 34) indicate that restraint against 
edge rotation and longitudinal movement significantly increases the buckling load, not 
enough is known about the edge restraint of actual cylinders to warrant taking 
advantage of these effects unless substantiated by tests. 

For layered or unstiffened orthotropic cylindrical shells, the available test data are 
quite meager (refs. 36 and 37). For configurations where the coupling coefficients C x , 
Cy, C X y, and K X y can be neglected, it is recommended that the buckling load be 
calculated from the equation 
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The correlation factor y is taken to be of the same form as the correlation factor for 
isotropic cylinders [eq. (5)] with the thickness replaced by the geometric mean of the 
radii of gyration for the axial and circumferential directions. Thus 

7 = 1 - 0.901 (1 - e'0) (45) 


where 
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( 46 ) 


4.3.2 Bending 


Theoretical and experimental results (refs. 24, 29, and 38 to 40) indicate that the 
critical maximum load per unit circumference of a stitfened cylinder in bending can 
exceed the critical unit load in axial compression. In the absence of an extensive 
investigation, it is recommended that the critical maximum load per unit circumference 
of a cylinder with closely spaced stiffeners be taken as equal to the critical load in axial 
compression, which is calculated from equation (37) multiplied by a factor of 0.75. 


For layered or unstiffened orthotropic cylinders with negligible coupling coefficients, 
it is recommended that the maximum unit load be calculated by means ot 
equation (44) with 7 replaced by 


7=1- 0.731 (1 -e'4>) 


where 



(47) 


(48) 


4.3.3 External Pressure 


The counterpart of equation (37) for stiffened orthotropic cylinders under lateral 
pressure is given by 



A 

‘ 1 1 

A 1 2 A 1 3 


A 21 

A 2 2 A 2 3 

r 

A.„ 

A 32 A 33 

7 2 


An 

Al2 




An 

A 22 



( 49 ) 


19 



For hydrostatic pressure, the quantity n 2 shown in equation (49) is replaced by 


In the case of lateral pressure, m is equal to unity while n must be varied to yield a 
minimum value of the critical pressure, but is not less than 2. In the case of hydrostatic 
pressure, the value of m should be varied as well as n. For long cylinders, equation (49) 
is replaced by 


P = 



(50) 


If the coupling coefficients C x , Cy, C X y, and K X y can be neglected, the critical 
buckling pressure can be approximated by (ref. 22): 
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Equation (49) has been investigated primarily for isotropic cylinders with ring 
stiffeners (refs. 41 to 43). For closely spaced ring stiffening, references 41 and 42 show 
that the effectiveness of inside or outside rings depends on the shell and ring 
geometries. Generally, for shells with values of Z less than 100, outside rings are more 
effective than inside rings, while for values of Z greater than 500, the reverse is true. As 
the ring geometry varies, the effectiveness of outside stiffening tends to increase as the 
stiffness of the rings relative to the shell increases. Somewhat lower buckling pressures 
are given by the extremely complex but more accurate theory of reference 44, but the 
differences are not so significant as to warrant its use. 

The experimental results for ring-stiffened cylinders described in reference 45 are in 
reasonably good agreement with the theoretical results of equation (49). For cylinders 
of all types, it is recommended that the buckling pressure calculated from 
equation (49) be multiplied by a factor of 0.75 for use in design, as has been 
recommended for unstiffened isotropic cylinders of moderate length. 
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4.3.4 Torsion 


The problem of torsional buckling of orthotropic cylinders has been treated in 
references 22 and 46, which do not take coupling between bending and extension into 
account, and in reference 47, which does. If coupling effects are negligible, the critical 
torque of moderately long cylinders can be estimated from the relationship (ref. 22): 
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Reference 47, however, shows that coupling effects are quite important for cylinders 
stiffened by closely spaced rings. For long shells, internal rings are generally more 
effective than outside rings; for short shells, the reverse is true. In the absence of 
general formulas or graphs to cover the entire range of parameters that should be 
considered, the equations of reference 47 should be solved for each specific case 
considered. 

The test data of reference 48 are in good agreement with theoretical predictions but 
are insufficient to provide an adequate test of the theory. It is therefore recommended 
that theoretical critical torques be multiplied by a factor of 0.67 for moderately long 
cylinders. 


4.3.5 Combined Bending and Axial Compression 

On the basis of theory (refs. 24, 29, and 38) and limited experimental data (refs. 28 
and 29), a straight-line interaction curve is recommended for orthotropic cylinders 
subjected to combined bending and axial compression. The critical combinations of 
loading are thus given by 


Rc + Rb = l 


(55) 
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4.3.6 Elastic Constants 


The values of the various elastic constants used in the theory of buckling of 
orthotropic cylinders are different for different types of construction. 


4.3. 6.1 Stiffened Multilayered Orthotropic Cylinders 

Some widely used expressions for this type of cylinder are: 
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where the subscript k refers to the material and geometry of the k™ layer of an 
N-layered shell. See figure 7. A proper choice of the reference surface can make at least 
one of the coupling coefficients vanish. For example, if A is taken as 
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the coefficient C X y vanishes and if 
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the coefficient K X y vanishes. 
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4.3.6. 2 Isotropic Cylinders with Stiffeners and Rings 


For a cylinder consisting of a stiffened single isotropic layer and for a reference surface 
at the center of the layer, equations (56) to (66) reduce to 
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4.3. 6.3 Ring-Stiffened Corrugated Cylinders 


The following formulas are commonly used to calculate the required stiffnesses of 
ring-stiffened corrugated cylinders, with the choice of formula depending on the 
different assumptions which may be made: 
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Slightly different stiffnesses are given in reference 39. 

4.3.6.4 Waffle-Stiffened Cylinders 

Stiffnesses for cylinders with waffle-like walls are described in references 49 to 5 1 . 

4.3.6 5 Special Considerations 


In some designs of stiffened cylinders, the skin may buckle before failure of the 
cylinder. Buckled sheet is less stiff than unbuckled sheet. The decreased stiffness can 
be calculated by methods similar to those presented in references 3 1 , 40, and 52. 
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4.4 Isotropic Sandwich Cylinders 


The term “isotropic sandwich” designates a layered construction formed by bonding 
two thin lacings to a thick core. Generally, the thin facings provide nearly all the 
bending rigidity of the construction. The core separates the facings and transmits shear 
so that the facings bend about a neutral axis. The core provides the shear rigidity of the 
sandwich construction. 

Sandwich construction should be checked for two possible modes of instability failure: 
(1) general instability failure where the shell fails with core and facings acting together, 
and (2) local instability failure taking the form of dimpling of the faces or wrinkling of 
the faces. See figure 8. 
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Figure 8 

Types of failure of sandwich shells 


If the isotropic sandwich shell has thin facings and the core has relatively little bending 
stiffness, then for unequal thickness facings, the bending stiffness is given by 
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and for equal thickness facings by 
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The extensional stiffness for unequal thickness facings is given by 
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The transverse shear stiffness for an isotropic core is given by 


and for equal thickness by 
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The stiffness of other types of sandwich construction are given in references 53, 54, 
and 55. 


4.4.1 Axial Compression 


Investigations of buckling behavior of isotropic sandwich circular cylinders in axial 
compression are reported in references 21 and 56. Design information from these 
references is given in figures 9 and 10. 

Figure 10 is the more convenient of the two figures to use; it is applicable to all but 
unusually short cylinders [7Z < jr 2 / (1 + R)]. Figures 9 and 10 are based on the 
small-deflection buckling theory and should be used in conjunction with the 
correlation factor of figure 1 1 to predict buckling loads. Figure 1 1 is based on 
equation (45), given earlier for orthotropic cylinders. For the present application, the 
parameter 0 becomes 



This procedure is consistent with the procedures given earlier for other types of 
construction when shearing of the core does not contribute significantly to the 
buckling deformations; that is, when N 0 /Dq of figure 10 is small. As shearing 
deformations become more pronounced, the correction resulting from application of 
the factor 7, as prescribed above, decreases and becomes zero in the limiting condition 
of buckling from a weak core [(N 0 /Dq) > 2] . 
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Buckling coefficient 



7Z 


Figure 9 

Buckling coefficients for simply supported isotropic sandwich 
circular cylinders subjected to axial compression G xz /G yz = 1.0 



N 0 /D q 

Figure 10 

Buckling of moderately long, simply supported, isotropic sand- 
wich circular cylinders subjected to axial compression 
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Correlation factor. 
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Figure 11 

Correlation factors for isotropic sandwich circular cylinders subjected to axial 
compression 


A weight-strength study based on figure 1 0 and published values for the shear stiffness 
of honeycomb cores (ref. 57) indicate that unusually lightweight cores are more 
desirable than heavier cores. Until adequate test data are obtained to substantiate this 
indication, however, designs should be limited to sandwiches with rather heavy cores 
(5 >0.03). Sandwich plates with light honeycomb cores are susceptible to additional 
modes of deformation, and failure may result from intracell buckling, face wrinkling, 
or an interaction of one or both of these modes with a cylinder-buckling mode. In 
addition, small buckle-like deformations have been known to occur in actual structures 
long before the theoretical buckling load is reached. (See, for example, ref. 19, p. 217.) 
This behavior requires that the structure be capable of resisting internal moments and 
shears in addition to the directly applied loads. In the case of sandwich cylinders, the 
moments and shears may cause core buckling or shear failure of the core. 


The only known method of preventing these core failures is to use relatively heavy 
cores which have considerable strength in crushing and shear. Some guidance as to how 
heavy the cores should be can perhaps be gleaned from the bending tests that have 
been made on multiweb beams. The internal structure of these beams is subject to the 
same types of loads, if not to the same degree, as the cores of loaded sandwich plates. 
Reference 58 indicates that honeycomb cores with a density ratio 5 = 0.03 should be 
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adequate to prevent core failure. Large margins against failure in intracell buckling and 
wrinkling can be obtained with rather heavy cores (5 <0.03) with little or no weight 
penalty. Moreover, when heavy cores are used, approximate equations are adequate for 
predicting failures in the intracell buckling and face wrinkling modes. The following 
equations may be used for this purpose. For intracell buckling (refs. 55 and 59): 
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where S is the core cell size expressed as the diameter of the largest inscribed circle and 
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where E and E tan are the elastic and tangent moduli of the face-sheet material. If 
initial dimpling is to be checked, the equation 
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should be used. The sandwich will still carry load if initial dimpling occurs. Critical 
wrinkling stresses are predicted by references 19 (p. 341) and 55. 


o x = 0.50 (E sec E z G xz ) - 1 (95) 

where E z is the modulus of the core in a direction perpendicular to the core and G xz is 
the shear modulus of the core in the x-z plane. If biaxial compressive stresses are 
applied to the sandwich, then the coefficients of the equations must be reduced by the 
factor ( 1 + f 3 )‘“ (ref. 60) where 

minimum principal compressive stress in lacings (96) 

* maximum principal compressive stress in lacings 


Wrinkling and intracell buckling equations which consider strength of bond, strength of 
foundation, and initial waviness of the facings are given in references 55,61, and 62. 

The plasticity correction factor given by equations (7) and (8) for isotropic cylinders in 
axial compression also may be applied to isotropic sandwich cylinders. The factor is 
applicable to sandwich cylinders with stiff cores and becomes somewhat conservative 
as the shear stiffness of the core is decreased (ref. 58). 
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4.4.2 Bending 


The buckling equations given in Section 4.2.1 for circular cylinders in axial 
compression may be used for cylinders in bending, provided the correlation factor 7 is 
taken from figure 1 2 instead of from figure 1 1 . Figure 1 2 is based on equation (47), 
given earlier for orthotropic cylinders in bending. 



Figure 12 

Correlation factors for isotropic sandwich circular cylinders 
subjected to bending 


4.4.3 Lateral Pressure 


A plot of ky against 7Z, constructed from the data given in reference 63, is given in 
figure 13. The straight-line portion of the curve of figure 13 for a sandwich cylinder 
with a rigid core (R = 0) is given by the equation 
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There are no experimental data to substantiate figure 13; experience with isotropic 
cylinders, however, suggests that a factor 7 equal to 0.56 is to be used with this figure. 
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Figure 13 

Buckling coefficients for simply supported isotropic sandwich 
circular cylinders subjected to lateral pressure G^/Gy; =1.0 

Here, as with sandwich cylinders in axial compression or bending, designs should be 
limited to sandwich cylinders for which the density ratio 5 is 0.03 or greater, unless 
the design is substantiated by adequate tests. 

The plasticity factors for isotropic cylinders subjected to external pressure, expressed 
by equations (22), (23), and (24), may be used for isotropic sandwich cylinders 
subjected to lateral pressure. 


4.4.4 Torsion 


Isotropic sandwich cylinders in torsion have not received the same attention as 
cylinders in compression. Rigid- and weak-core cases are reasonably well defined. While 
the transition region between rigid and weak cores is not as well defined, it is probably 
sufficient for design purposes. Information on the transition region is given in 
references 53 and 63, the latter of which was used to construct the plot of figure 14, 
which applies to sandwich cylinders with cores exhibiting isotropic shear behavior 
G xz /Gy Z = 1 . The curves of this figure are discontinuous at the value of 7 Z where the 
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Figure 14 

Buckling coefficients for simply supported isotropic sandwich 
circular cylinders subjected to torsion G xz /Gy Z = 1.0 

buckling coefficient k xy becomes equal to 1/R, indicating a change of mode of 
buckling at that point. Reference 53 does not support this behavior, but it does not 
cover a sufficiently wide range of geometric proportions to be used in the construction 
of the figure. In addition, reference 53 indicates that there was some scatter in the 
calculated results used to construct the charts of that reference. In the ranges where 
comparisons between the data of references 53 and 63 could be made, only rather 
small discrepancies were noted. The straight-line portion (yZ > 170) of the curve of 
figure 14 for a rigid core (R = 0) is given by the equation 

N XV C 2 4 

k xv = — = 0 34 (tZ ) 4 (98 > 

' 3 7T“ L)| 

Experimental data are not available to substantiate figure 14 for most sandwich 
cylinders. Experience with isotropic cylinders suggests that 0.586 is the factor y to be 
used with the figure. Here, as with sandwich cylinders for which the density ratio of 5 
is 0.03 or greater, the same factor should be used unless the design is substantiated by 
adequate tests. 

The plasticity factor given by equation (29) for isotropic cylinders in torsion may be 
used also for isotropic sandwich cylinders in torsion. 
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4.5 Cylinders with an Elastic Core 


The term “cylinder with an elastic core” defines a thin cylindrical shell enclosing an 
elastic material that either can be solid or have a hole in its center. This type of shell 
closely approximates a propellant-filled missile structure. The propellant is generally of 
a viscoelastic material and therefore is strain-rate sensitive. The core modulus should be 
obtained from tension or compression tests of the core material simulating its expected 
strain rate. 

Although there are some analytical data for orthotropic shells (ref. 64), design curves 
are given only for isotropic shells and cores. The inverse problem of a core or cushion 
on the outside of the cylindrical shell is analyzed in reference 65. Not enough data are 
available, however, to recommend design curves for this problem. 


4.5.1 Axial Compression 


The buckling behavior of cylindrical shells with a solid elastic core in axial compression 
is given in reference 66. Analytical results obtained from this reference are shown 
graphically in figure 15. For small values of < 
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This approximatien is accurate for 0j less than 1/2. For larger values of <f> j , say <p 1 
greater than 3, 
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The experimental data tabulated in reference 66 suggest that the value of 7 to be used 
in calculating o c can be taken as that for isotropic cylinders in compression. Then 

7=1- 0401 (I — e"0) (103) 
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Figure 15 

Variation of compressive buckling stress with core stiffness parameter 


where 

♦-T sTi (,04) 

The plasticity correlation factors given by equations (7) and (8) for isotropic cylinders 
in axial compression may be applied also to the cylinder with an elastic core. 


4.5.2 External Pressure 

Analytical curves for the lateral pressure case are presented in reference 66. A plot of 
kp C against ?j- for -I = 100, 200, 500, or 1000 is shown graphically in figure 16. The 
parameter kp C is expressed by 

kpc = £ (105) 

These curves are to be used for finite cylinders loaded by lateral pressures. 
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Figure 16 

Variation of buckling pressure coefficient with length and modulus 
ratio (/x= 0.3, n c = 0.5) 
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Some cylinders are long enough for the critical pressure to be independent of length 
(fig. 16); the single curve shown in figure 17 can then be used. The straight-line portion 
of the curve can be approximated by the equation 
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The few experimental data points available indicate good agreement between analysis 
and experiment, but one test point falls 4% below theory. Hence, a correlation factor 
of 0.90 is recommended for use in conjunction with the curves in figures 16 and 17. A 
reinvestigation of this factor may be warranted as more data become available. 


The plasticity factors given by equations (22) to (24) for isotropic cylinders subjected 
to external pressure may be used also for cylinders filled with an elastic core. 
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Figure 17 

Buckling pressure coefficients for long cylinder with 
a solid core 
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4.5.3 Torsion 


The buckling behavior of cylindrical shells with an elastic core is analytically described 
in reference 67 and is shown graphically in figure 18. 

For small values of 0 3 ,0 3 <7, the analytical results can be approximated by 


— = 1+ 0.1603 (108) 
r cr 


where 



(109) 


and r cr is the torsional buckling stress given by equation (26), with 7 equal to unity. 
When 0 3 is greater than 10, the analytical results follow the curve 


r 2 

— ■ = 1 + 0.25 03 4 

T er 


( 110 ) 



Figure 18 

Torsional buckling coefficients for cylinders with an elastic core 
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Experimental data are not available for this loading condition. The experimental points 
obtained for cylinders with an elastic core for axial compression and external pressure, 
however, show better correlation with theory than the corresponding experimental 
results for the unfilled cylinder. Hence, conservative design curves can be obtained by 
calculating T cr in equations (108) and (110) with the correlation and plasticity factors 
for isotropic cylinders of equations (28) and (29). 


4.5.4 Combined Axial Compression and Lateral Pressure 


Interaction curves for cylinders with an elastic core subjected to combined axial 
compression and lateral pressure are shown in figure 19. These curves were obtained 
analytically in reference 66 and indicate that for a sufficiently stiff core, the critical 
axial compressive status is insensitive to lateral pressure and, similarly, the critical 
lateral pressure is insensitive to axial compression. Until more experimental data 
become available, the use of a straight-line interaction curve is recommended for 
conservative design. 



Figure 19 

Interaction curves for cylinders (r/t = 300) with an elastic core 
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4.6 Design of Rings 


Little information is available on which to base the design of rings for cylinders not 
intended to be subject to general instability failures. The criterion of reference 68 is 
frequently cited as applicable to cylinders subjected to bending or compression. 
Unfortunately, this criterion is empirical and is based on data from test cylinders with 
proportions of little interest in contemporary design. A few checks made on cylinders 
in use have indicated that the criterion usually is conservative, but this may not be so 
in certain cases. See references 29 and 69. 

A less direct procedure for designing rings may be used. It consists simply of 
calculating the failing load of the cylinder in the so-called general-instability mode, 
which involves failure of the rings, as well as calculation of the failing load of the 
cylinder for wall failure between rings. Both calculations are made for several ring 
weights. If such calculations are plotted against ring weight, the weight necessary to 
force failure in the desired mode can be ascertained. In addition, the amount of error 
in weight from uncertainties in the calculations can be judged. Presumably, there may 
be some interaction between failing modes; thus, somewhat heavier rings than those 
indicated by the calculations should be used. 

This method of designing rings is, of course, applicable to all types of loading and to all 
types of wall construction. It also has the advantage of giving the designer some feeling 
for the influence of the various factors which determine ring weight. 

A study of references 69 and 70, which present general linear analyses of ring-stiffened 
isotropic cylinders in torsion and of orthotropic cylinders in compression, indicates 
that the recommended procedure gives the same result as general theory for all 
cylinders except those with a single ring dividing the cylinder into two equal bays. 
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SP-8001 

(Structures) 

Buffeting During Launch and Exit, May 1964 

SP-8002 

(Structures) 

Flight-Loads Measurements During Launch and Exit, 
December 1964 

SP-8003 

(Structures) 

Flutter, Buzz, and Divergence, July 1964 

SP-8004 

(Structures) 

Panel Flutter, May 1965 

SP-8005 

(Environment) 

Solar Electromagnetic Radiation, June 1965 

SP-8006 

(Structures) 

Local Steady Aerodynamic Loads During Launch and 
Exit, May 1 965 

SP-8008 

(Structures) 

Prelaunch Ground Wind Loads, November 1965 

SP-8009 

(Structures) 

Propellant Slosh Loads, August 1968 

SP-8010 

(Environment) 

Models of Mars Atmosphere (1967), May 1968 

SP-8011 

(Environment) 

Models of Venus Atmosphere (1968), December 1968 

SP-8012 

(Structures) 

Natural Vibration Modal Analysis, September 1968 

SP-8014 

(Structures) 

Entry Thermal Protection, August 1968 

SP-8015 

(Guidance and 
control) 

Guidance and Navigation for Entry Vehicles, 
November 1968 
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